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(57) Abstract: 

PROBLEM. JTO BE SOLVED: To provide integration 
system/method effective even for highly dimensional 
Monte Carlo calculation, and high in convergence speed 
without impairing calculation time. 

SOLUTION: One-dimensional quantity is calculated from 
multidimensional random numbers. It is desirable that 
one-dimensional quantity is adjusted to the feature of a 
function to be integrated. One-dimensional quantity is 
classified into several groups by a value and the number 
of generation times is recorded at every group. The 
accumulation value of the value of the function to be 
integrated, which corresponds to the multidimensional 
random number at that time, is recorded in accordance 
with previous grouping. The operation is repeated for N 
times and a ratio between the distribution of 
one-dimensional quantity, which is logically predicted 
for the respective groups, and the recorded number of 
generation times is recorded as weight for the group. 
The product of the accumulation value of the integrated 
function and weight is calculated at every group and the 
sum of all the groups is taken. The sum is divided by 
the sum of weight against all the groups. Then, an 
integration value by corrected Monte Carlo analysis is 
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void gettheory (double theory [STATNUM] ) ; 

double getaverage (double theory [STATNUM] , 

int ustat [STATNUM], double p sum [STATNUM], long cnt) ; 

int 

main 0 
{ 



long 
double 
double 
int 
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u[D]; 

norm, price, priceave; 
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double psum [STATNUM], theory [STATNUM] ; 

for (i = 0; i < STATNUM; i++) 
{ 

ustat[i] = 0; 
psum[i] = 0; 

} 

gettheory (theory) ; 

for (jj = 1; jj <= maxcount; jj++) 

{ 

getrandom(u) ; /* get D-dim random sequences */ 

norm = calcnorm(u); /* calc norm of random sequence */ 
nindex = (int) ( (norm+STATNUM*STATLEN / 2) / STATLEN) : 
price = value (u); /* get a price from a senario */ 
ustat [nindex] ++ ; 
psum[nindex] += price; 

} 

priceave = getaverage (theory, ustat, psum, maxcount); 
printf ("result : %. 16f¥n", priceave) ; 

} 

void gettheory (double theory [STATNUM] ) 
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int mid = STATNUM / 2; 
int i, j; 

double dx = STATLEN / TNUM; 
double c, f, sum, x, total = 0; 

c = dx / sqrt(2 * D * 3. 1415926536); 

for (i = 6; i < mid - 1; i++) 

{ 

sum = 0; 

for (j = 0; j < TNUM; j++) 
{ 

x = i * STATLEN + ( j + 0. 5) * dx; 
f = exp(- x * x / (2 # D)) * C ; 
sum +- f; 

} 

theory [mid + i] = theory [mid - i -1] = sum; 
total += sum; 

} 

theory [0] = theory [STATNUM -1] = 0. 5 - total; 

} 

double getaverage (double theory [STATNUM] , 

int ustat [STATNUM], double psum [STATNUM], long cnt) 

{ 

double pricesum = 0.0, wsum = 0.0, weight; 
int i; 

for (i = 0; i < STATNUM; i++) 
{ 

weight = theory [i] / ustat [i] ; 
pricesum += psum[i] * weight; 
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return pricesum / wsum; 
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